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We analyze the smallest Dirac eigenvalues by formulating an effective theory for 
the QCD Dirac spectrum. We find that in a domain where the kinetic term of the 
effective theory can be ignored, the Dirac eigenvalues are distributed according to a 
Random Matrix Theory with the global symmetries of the QCD partition function. 
The kinetic term provides information on the slope of the average spectral density 
of the Dirac operator. In the second half of this lecture we interpret quenched 
QCD Dirac spectra at nonzero chemical potential (with eigenvalues scattered in 
the complex plane) in terms of an effective low energy theory. 

1 Introduction 

There are two essentially different approaches to many-body problems. The first 
approach is to study them by means of Monte-Carlo simulations of the microscopic 
theory. The second approach is to isolate the relevant degrees of freedom and 
to describe them by an effective theory. Both approaches have their merits and 
generally their complementarity leads to a deeper understanding of the underlying 
phenomena. In this lecture we will focus on Quantum Chromo-Dynamics (QCD) 
which is the theory of the strong interactions. A.-great deal of effort has been 
devoted to Monte Carlo simulations of lattice QCD a. They provide a firm footing 
for our understanding of nonperturbative phenomena such as confinement and chiral 
symmetry breaking. For QCD at low energy an alternative approach is possible. 
Because of spontaneous breaking of chiral symmetry, QCD at low energy reduces to 
a theory of weakly interacting Goldstone bosons. Although this theory cannot be 
derived from QCD by means of an ab-initio calculation, its Lagrangian is determined 
uniquely by chiral symmetry and Lorentz invariance 0. The validity of this low- 
energy theory is based on the presence of a mass-gap which is a highly nontrivial 
and nonperturbative feature of QCD. 

Chiral symmetry is spontaneously broken at low temperatures. Mainly through 
lattice simulations, it is now widely accepted that a chiral restoration transition 
takes place at a temperature of about 140 MeV. The order parameter of this tran- 
sition is the density per unit volume of Dirac eigenvalues near zero. This is the 
reason why we pare interested in the properties of the smallest eigenvalues of the 
Dirac operator era. The situation at finite baryon density and zero temperature is 
much less well understood. Monte-Carlo simulations are only possible if the fermion 
determinant is ignored. It has been shown that the so C^Ued quenched approxima- 
tion fails spectacularly at nonzero chemical potential an. On general grounds it 
is expected that a chiral phase transition occurs at a value of the chemical poten- 
tial where it becomes advantageous to create the lightest particles with nonzero 
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baryon number. In QCD these are the nucleons, but in the quenched approxima- 
tion, they are Goldstone modes made from quarks and conjugate anti-quarks a. At 
nonzero chemical potential, the relation between the chiral order parameter and 
the QCD Dirac spectrum, which is now scattered in the complex plane, is much 
less transparent. The failure of quenching can be understood as the absence of 
spectral ergodicity; at finite density the ensemble averaged Dirac spectrum is com- 
pletely different form the spectral averaged Dirac spectrum. We will investigate 
the Dirac spectrum by means of an effective theory. Our approach is similar to 
the one used for QCD with two colors and for QCD with adjoint quarks Era. By 
investigating the properties of the smallest Dirac eigenvalues we hope to obtain a 
better understanding of this problem and other related issues. 

In this lecture we wish to discuss to what extent spectra of the QCD Dirac 
operator can be derived from an effective low energy partition function. In the 
first half of this talk we formulate an effective theory for the Dirac spectrum based 
on a spectrum generating function that in addition to the usual quarks contains 
bosonic ghost quarks to properly normalize the spectral density. This trick has. 
widely used in the super-symmetric formulation of Random Matrix Theory BE 
We identify a domain where the results for the resolvent are given by a chiral 
Random Matrix Theory with the symmetries of the QCD partition function. We 
find that our results are in agreement with recent lattice simulations. In the second 
half of this lecture we discuss quenched QCD Dirac spectra at nonzero chemical 
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2 The QCD partition function 

The QCD partition function for Nf quarks with mass m, temperature T and quark 
chemical potential \i is given by 

h qcd-c n 

Z(m,fx,T) = Tre * , (1) 

where Hqcv is the Hamiltonian of QCD and N is the quark number operator. The 
trace is over all states of the theory. This partition function can be rewritten as a 
Euclidean Feynman path integral 

Z(m, (j,, T) = (det N f (D + m)) YM , (2) 

where the average of the fermion determinant is over the Euclidean Yang-Mills 
action and the Dirac operator is denoted by D. In this lecture we focus on the 
basic symmetry properties of D and its explicit representation is not necessary. For 
simplicity we take all Nt quark masses equal to m. 

2. 1 The eigenvalues of the Dirac operator 

The main subject of this talk are the properties of the smallest eigenvalues of 
the QCD Dirac operator. In a chiral representation of the gamma-matrices, the 
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Euclidean Dirac operator has the block structure 

where d is the covariant derivative of the color group. For \i ^ the Dirac operator 
does not have any hcrmiticity properties and the eigenvalues will be scattered in 
the complex plane. We will return to this case in the second half of this talk. For 
/i = the Dirac operator is anti-Hermitian, = — D and has purely imaginary 
eigenvalues. The nonzero eigenvalues occur in pairs HXk, whereas zero eigenvalues 
related to the topology of the gauge fields remain unpaired. For simplicity we will 
restrict the discussion in this lecture to gauge fields with a trivial topology where 
all eigenvalues occur as pairs. For example, our gauge fields cauld be composed 
of fields of an equal number of instantons and anti-instantons □. In terms of its 
eigenvalues, the Euclidean partition function can be written as 

Z(m )t i ) T) = (T[(i\ k + m) N n . (4) 

\ k I YM 

The order parameter of the chiral phase transition is the chiral condensate 
£ = lim lim -tt^-Sto log Z 

m^O V^od VNf 

= lim lim 77^/ — — ) = lim lim i / — 2w „ \ ,(5) 

m^o ^00 V Y \ m + i\ k / QCD m^o y^oo V ^ Q \ m 2 + \\ / QCD 

where the average with label QCD includes both the fermion determinant and the 
Yang-Mills action. At finite volume, the chiral condensate is zero for m = 0. Only 
if the thermodynamic limit is taken before the chiral limit (m — » 0) can the chiral 
condensate become nonzero. If this is the case, chiral symmetry is broken sponta- 
neously. The situation is analogous to the magnetization in a Heisenberg model. At 
finite volume and zero magnetic field the magnetization is zero. Spontaneous mag- 
netization arises if the thermodynamic limit is taken before putting the external 
field to zero. According to Goldstone's theorem spontaneous breaking of a con- 
tinuous symmetry implies the existence of massless Goldstone bosons. Because of 
confinement QCD has a mass gap and, at low energy, the partition function is dom- 
inated by the Goldstone modes. The effective theory describing their interactions 
then follows from chiral symmetry and Lorentz invariance. 

One of the questions we would like to address in this lecture is whether there is a 
relation between the existence of a well-defined low-energy theory and the spectrum 
of the Dirac operator. A second, seemingly unrelated question is the connection of 
this low-energy with Random Matrix Theory. 

3 Symmetries of the Spectrum Generating Function at fi = 

We will study the Dirac spectrum by means of the resolvent defined by 

T,( z ) = -(Ti— — ) . (6) 

V) V\ D + z/ QCD [) 
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The mass z in the resolvent is not related to the mass m in the fermion determinant. 
Such mass is known in the literature as a valence quark mass. However, we will use 
the more appropriate name of spectral mass. 
In terms of the spectral density 



k 

the resolvent can be written as 

m = I / Pi x )d x-± 



+ z 



(7) 



(8) 



Therefore, ~S(z) is an analytic function with a cut on the imaginary axis. This 
identity can be inverted by taking the discontinuity across the cut 



^l = lim^(E(iA + e)-S(iA-e)). 



(9) 



In order to obtain a generating function for Il£zLwe use a method that has been 
widely used in the theory of disordered systems laO, namely ElHl§tL3 



Z sp (z,J,m) = ( dct } D + Z+ ; J) dct N f(D + m)) YM . 



det(D + z) 



The resolvent is then given by 



V 



djZ(z, J, m) 



(10) 



(11) 



,7=0 



The spectrum generating function (10) contains Nf + 1 fermionic quarks and one 
bosonic quark. In addition to the chiral symmetry, this partition function also 
contains a super-symmetry that mixes fermionic and bosonic quarks. This can be 
seen by rearranging the fermion determinant as 



det 



m id 



det 



id m 



We observe that for m — z 
(id \ 



id 



i<# 



V 



id) m J ym id^ 

J = the partition function is invariant under 
(id \ 



(12) 



U 



V 



id^ J 



id 



iS 



\ 



u- 



y- 



(13) 



id) J 



where U an V are (Nf + 1|1) x (Nf + 1|1) super-matrices. Mathematically, this 
symmetry group is known as Gl(Nf + 1|1) x Gl(Nf + 1|1). A G((l) subgroup is 
broken by the anomaly. The chiral condensate S = (ipip) is only invariant under the 
diagonal subgroup with U = V. Therefore the symmetry is broken spontaneously 
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to Gl(Nf + 1|1). As is the case in QCD, the masses of the Goldstone modes are 
given by the Gell-Mann-Oakes-Renner relation 

\/2m£ y/(m + z)E a/2zE 

where F is the pion decay constant. The Goldstone modes corresponding to a 
fermionic and a bosonic quark are fermionic whereas all other Goldstone modes are 
bosonic. 



4 Effective Low-Energy Theory 

The manifold Gl(Nf + 1|1) is not Riemannian, and is therefore not suitable as a 
Goldstone manifold. The Goldstone manifold is given by the maximum Riemannian 
submanifold of the symmetric superspace Gl(Nf + 1|1) which will be denoted by 
Gl(N f + 1|1) @0fl If we ignore certain complications related to the topological 
structure of the QCD vacuum, this partition function is given by 

Z{m,J,z) = f dUe -J^[^stra^u-^^ S t r M(u + u-^^ (15) 

JGl(N+f+l\l) 

where the mass. matrix M = diag(m, • • • , to, z + J, z).The super- matrix U is param- 
eterized as Ej'Ej 

U=( V p °)=eW. (16) 

Here, V is a U(Nf + l)-matrix, a and (3 are Grassmann valued vectors of length 
Nf + 1 and s is a real number. 

In order to estimate the relative importance of the two terms in the effective 
Lagrangian we expand the fields to second order in the pion fields $ = Tr a t a , 

C eS = \d^d^« + \Ml*l (17) 

where M a is one of the Goldstone masses given in (|l4|). In a box of volume L 4 , 
the smallest nonzero momenta are of the order ~ 1/L. Let us consider QCD 
in the chiral limit (with to = 0). Then the regular mesons are massless whereas 
mesons containing one or two spectral quarks have a nonzero mass given by (|l4|). 
Therefore, if JB 

zT. 1 

F^«J2 ( 18 ) 

the correlation functions involving spectral quarks are dominated by contributions 
from the zero momentum Goldstone modes. Nonperturbatively, the partition func- 
tion reduces, to a, group integral in this limit. This integral has been calculated 
analytically ll3'E3 resulting in the following dependence of the condensate on the 
spectral mass 

S(z) 
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As already could be observed from the generating function, it depends on z only 
through the combination fj, z = zVYj. This result was first obtained 113 by means of 
chiral Random Matrix Theory to be discussed in the next section. 

QCD is not the only theory that reduces to this effective partition function. In 
fact, any partition function with the same pattern of chiral symmetry breaking and 
a mass gap reduces to the_same low-energy theory. Explicit examples are given 
by the random flux-model Ej and other models with a hopping term, disorder and 
chiral symmetry I25EJ. A natural question is what is the simplest theory with the 
same zero momentum sector as QCD. This theory is chiral Random Matrix Theory. 



5 Chiral Random Matrix Theory 

In the sector of topological charge v and for Nf quarks with mass m, the chiral 
random matrix partitiou_function with the global symmetries of the QCD partition 
function is defined by qEj 

W =/wndet(; t ^)e-^^, (20) 

where W is anx (n+|^|) matrix and N = 2n+\v\. As is the case in QCD, we assume 
that the equivalent of the topological charge v does not exceed \[N , so that, to a 
good approximation, n — N/2. Then the parameter E can be identified as the chiral 
condensate and N as the dimensionless volume of space time (Our units are defined 
such that the density of the modes N/V =1). The matrix elements of W are either 
real ((3 = 1, chiral Gaussian Orthogonal Ensemble (chGOE)), complex (f3 = 2, 
chiral Gaussian Unitary Ensemble (chGUE)), or quaternion real ((3 = 4, chiral 
Gaussian Symplectic Ensemble (chGSE)). For QCD with three or more colors and 
quarks in the fundamental representation the matrix elements of the Dirac operator 
are complex and we have (3 = 2. The ensembles with (3=1 and [3 = 4 are relevant 
in the case of two colors and adjoint fcrmions, respectively. The reason for choosing 
a Gaussian distribution of the matrix elements is its mathematical simplicity. It can 
be shown that the correlations of the eigenvalues on the scale-oij-the Jey£l snacing-do 
not depend on the details of the probability distribution BNkaMyHayMEa 
Together with_the Wigner-Dyson ensembles and the superconducting random 
matrix ensembles till the chiral ensembles-can be classified according to the Cartan 
classification of large symmetric spaces E3. 



6 Scales in the Dirac Spectrum 

For a nonzero value of the chiral condensate £ we can identify three important scales 
in the Dirac spectrum. The first scale is the smallest nonzero eigenvalue of the Dirac 
operator given by A min = l//o(0) = tt/Y>V. The second scale is the spectral mass 
for which the Compton wavelength of the associated Goldsto pe . bps ttas is equal to 
the size of the box. As discussed above this scale is given by BEaOEJ 

F 2 

E ° = m- (21) 
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In mesoscopic physics this scale is known as the Thouless energy. It is given 
by the inverse diffusion time of an electron through a sample of length L. In 
Euclidean QCD, this scale can be interpreted in terms of diffusive motion of quark 
in 4 Euclidean dimensions and one artificial time dimension Ej. A third scale is given 
by a typical hadronic mass scale. The three scales are ordered as A m i n C B c « A. 

For spectral masses z <C E c the kinetic term in the effective action can be ne- 
glected and the low-energy partition function can be reduced to a zero dimensional 
integral. As already mentioned, the results .fox £(z) obtained from the spectral 
partition function and from chRMT coincide lialla in this domain. 
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Figure 1. The spectral mass dependence of the chiral condensate T,(z) plotted as £(z)/£ versus 
zVT,. The dots and squares represent lattice results by the Columbia group EJ for values of /3 as 
indicated in the label of the figure. 




In the domain E c <C z <C A the kinetic term has to be taken into account. The 
slope of the Dirac spectrum at A = can be obtained from a one-loop calculation. 
For Nf massless flavors the result is given by eIEjO 

P '(0) _ (N f -2)(N f +f3) S 



p(0) 16tt/3 F 4 ' 



(22) 



Here, f3 denotes the Dyson index of the Dirac operator defined before. 

The domain below E c has been investigated extensively by means of lat- 

J^B^S^ T ltS ^ b6en 

found li^oaac3'C3Lllcacj' EiJEiltiToOt^l^tJea. In Fig. 1 we show a compar- 
ison of_the ratio E(^)/E versus zVT, of lattice results obtained by the Columbia 
group c3 and eq. Jl9| ) for Nf = and 2. The point where the lattice data depart 



from the chRMT result roughly coincides with the scale E c defined in eq. (21) 
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7 QCD Dirac spectra at /i ^ 



The QCD partition function at /i ^ is given by 



(23) 



Below we will derive the Dirac spectrum from the property that for T — > only 
states with E a /N a < (i contribute to the partition function. 

For (i^O the eigenvalues are scattered in the complex plane. The spectral 
density is then given by Ej 



P {z) = -d z .G(z), 
with the resolvent G(z) defined as usual by 

G(z) 



1 1 
— Tr 

V D + z 



det N f (D + ■ 



(24) 



(25) 



YM 



Writing out explicitly the real and imaginary parts, we observe that G(z) is the 
electric field in the plane of charges located at the positions of the eigenvalues. 
The spectrum of the Dira^^pep-tar-bas. Jaeee studied, .extensively in the context of 
Random Matrix Theory 163000^^^^ ■ |^6£j^,^i[ Below, we will obtain the 
quenched Dirac spectrum using an effective low energy. To lowest order in fi 2 our 
results agree with quenched Random Matrix Theory. 

At [i a new ingredient enters in the definition of the the spectrum gener- 
ating function. In order to assure convergence of the bosonic integral we need an 
additional factor det(Z?^ + z*) in the denominator, and a corresponding factor in 
the numerator, 



Z(z, m) 



I det(L> + z + J) det(£>t + z* + J*)det Nf (D + m) 
\ det(D + z)det(L>t + z*) 



(26) 



YM 



The additional factor can be interpreted in terms of conjugate quarks with a baryon 
number opposite to that of ordinary quarks This opens the possibility of baryonic 
Goldstone modes consisting of a quark and a conjugate anti-quark. For simplicity 
we only discuss the quenched case where Nf = 0. 

In order to construct the effective partition function we have to identify the 
symmetries of the spectrum generating function. To this end we rewrite the product 
of the determinant and the conjugate determinant as follows 



det(D + z)det(L> + z*) = dct 



det 



id + fi z 
z iS + 



det 



fid- 



\ 



[i z 
id — fi 
id^ 4 
-z* 



— id - 

z* 



— z* 
[i 
id) - 



z 

-iS 



(27) 



As can be observed from the sign of the chemical potential, the baryon number of 
the conjugate quarks is opposite to that of the usual quarks. 
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In this case the chiral symmetry is broken spontaneously according to 

Gl(2\2) x Gl(2\2) -> Gl(2\2). (28) 

This symmetry is broken explicitly by the term proportional to /i as 

Gi(l,l|l,l)xG/(l,l|l,l). (29) 

This gives rise to an additional term in the effective action. In order to maintain 
invariance under both the left and right symmetry groups we need at least a term 
of second order in fi and U. One easily verifies that to this order, the only static 
term with the correct symmetry properties is given by 

J d i xfj, 2 StrUT 3 U- 1 T 3 , (30) 

where the symmetry breaking matrix T3 is a diagonal matrix with nonzero matrix 
elements equal to 1 except on the positions corresponding to the conjugate quarks, 
where they are -1. For example, for Nj = we have 



n = 



f 1 \ 

-1 

1 

V -v 



(31) 



If we ignore complications related to the topology of the gauge field configurations, 
the static part of the effective partition function is given by 

p ±T;VStr MiU+U- 1 )- ±V fj, 2 F 2 Str UTsU-iTa ^2) 

U£Gl(2\2) 

where the mass matrix is given by M — diag(z + J, z* + J*, z, z*). The integral is 
over the maximum Riemannian submanifold of Gl (2\2). The relative coefficient of 
the two symmetry breaking terms is determined by the condition that the partition 
function should have a singularity when 2/i becomes equal to the mass of the lightest 
meson. The coefficient of the term ~ /i 2 has to be chosen such that the effective 
meson mass vanishes at 

guaranteeing a nonanlyticity of the /i-dependence at this point. The value of this 
coefficient can be obtained more elegantly by means of a gauge principle. This 
construction haa^been performed for QCD with two colors and for QCD with ad- 
joint fermions 0i3. In addition to the term proportional to /i 2 we then obtain the 
coefficients of the terms in the effective Lagrangian that are linear in fj, and the 
momentum. 

Baryonic Goldstone modes contain one ordinary quark and one conjugate quark 
both with mass z resulting in a square mass of 2RezE/F 2 . For /j 2 < RezE/2F 2 
only the vacuum state contributes to the QCD partition function, and thus 

Z(z,J) = e 2RcJSV \ (34) 
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In the effective theory this corresponds to the saddle-point U = 1. For /i 2 > 
Re zH/2F 2 , baryonic Goldstone modes contribute to the partition function result- 
ing in a nonzero baryon density. The density stays finite because of the repulsive 
interaction between the Goldstone modes. In terms of the effective partition func- 
tion this happens by the rotation of the saddle point from U = 1 to a nontrivial 
value. By making the Ansatz for the boson-boson and the fermion-fermion blocks 

Ubb = U F f = [ C ° S ® a Sm a ) , (35) 
V — sm cos 1 

one finds that cos — ERez /2fi 2 F 2 . A similar rotation of the saddle, point has 
been found in the analysis of nonhermitian random matrix models Ej^Bej. We then 
find the following J-dependence of the partition function 

„, . r VT, 2 ReJRez, 

Z(z, J) = exp[ ] • (36) 

The resolvent in both domains follows by differentiation with respect to the source 
J. For Rez > 2y 2 F 2 /Y 1 we find from (|J) 

G{z) = E and p(z) = 0, (37) 

and for Re z < 2/i 2 i* l2 /E the result for the resolvent following from ( |36| ) is given by 

E 2 Rez E 2 ,„_* 

We observe that the resolvent is continuous at the transition point. The eigenvalues 
are thus located inside a strip of width 4i^ 2 /K 2 /E. In terms of the interpretation 
of the resolvent as an electric field, an eigenvalue density that does not depend on 
the imaginary part of z results in a constant electric field outside of the strip of 
eigenvalues. Ignoring finite size effects this is indeed what is observed in quenched 
lattice QCD simulations 13. In the quenched approximation these results are in 
agreement with an explicit Random Matrix Model calculation. Therefore, to order 
fj, 2 the Random Matrix partition function reproduces exact QCD results. 

8 Conclusions 

We have shown that the correlations of QCD Dirac eigenvalues at a scale well 
below the hadronic mass scale can be obtained analytically. They are given by 
an effective theory for the generating function of the resolvent for the QCD Dirac 
spectrum. For mass scales below F 2 /EL 2 the kinetic term in the effective theory 
can be ignored resulting in eigenvalue correlations given by chRMT. Our results 
have been confirmed by numerous lattice QCD simulations. The same procedure 
can be followed for QCD at nonzero chemical potential. Our results indicate that in 
the quenched limit the eigenvalues are scattered in a strip with a width determined 
by the lightest meson of our theory. This shows that, to lowest nontrivial order in 
the chemical potential, chiral Random Matrix Theory provides an exact description 
of the QCD Dirac spectrum. 
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